Level correlations in disordered superconducting grains 



O 

o 
o 

(N 
0\ 



■4— » 
•4-* 



O 

o 



> 

O 
O 
O 

-i— > 
c3 



c 

o 
o 



X 



Revaz Ramazashvili 
Loomis Laboratory, University of Illinois at Urbana- Champaign, 
Urbana, IL 61801-3080, USA. 

I study the quasiparticle level correlations in a grain of a weakly disordered d-wave superconductor, 
and show that, in a wide intermediate energy range, they are characterized by a novel type of 
universal behavior. 
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Level correlations represent a fundamental property 
of an electron system. They reflect sensitivity of the 
quasiparticle spectrum to disorder and boundary con- 
ditions, and distinguish between extended, apd localized 
states-. Level correlations in metal grains,EJ~l!l mesoscopic 
wiresjj quantum dote!! SNS junctionsjj and supercon- 
ducting vortex coresETM have recently been a subject of 
much active study. Most simply, level correlations are 
quantified by the mean square deviation (SN E ) of the 
number Ne of levels in an energy interval of width E: 



(SN 



(Ne)Y 



where the angular brackets denote disorder average. 

In disordered systems, the energy levels of localized 
states are uncorrelated, and (SN E ) ~ (Ne^qBy con- 
trast, extended states are strongly correlated, 00 which is 
expressed in (SN E ) scaling as only a logarithm of (Ne)' 



(6N 2 E ) 



2C 



\n(N E ) < (N E ), 
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where 11 is the number of the quasiparticle diffusion 
modes.E3 Formula (0) not only measures the level cor- 
relations, but also shows that, in a disordered system, 
the level number variance is driven by diffusion modes. 
Formula ([!]) holds for E smaller than the Thouless en- 
ergy E c = D/L 2 , where D is the diffusion constant, and 
L is the system size. At these energies, the level cor- 
relations are remarkably universal, with the constant C 
being defined solely by the fundamental symmetries of 
the system. The low energy universality classes of disor- 
dered systems tfj^ been classified based on the symme- 
try arguments JI34L-3 

In this article, I study quasiparticle level correlations 
in a disordered grain of a d-wave superconductor in the 
presence of both the time reversal (T) and the spin ro- 
tation invariance (S). I show that, in a wide intermedi- 
ate energy range, the level correlations have the univer- 
sal form (Q), yet are different from those in a grain of 
a metal or of an s-wave superconductor (also invariant 
under T and S) - as well as different from lie level cor- 
relations in any of the previously chartecbj'Ej low energy 
universality classes. The main result is encapsulated in 



the number C of the diffusion modes. In a grain of nor- 
mal metal (or an s-wave superconductor) with both the 
time reversal and the spin rotation symmetries present, 
C = 4, since both the quasiparticle charge and the three 
components of the quasiparticle spin are conserved and 
propagate diffusively. By contrast, I show that, in a d- 
wave superconducting grain, the charge diffusion mode 
disappears, which leads to C — 3, despite the very same 
set of fundamental symmetries. This is a novel universal 
type of level correlations. 

The reason behind this result is that the impurity scat- 
tering leads to the exchange of charge between the quasi- 
particle subsystem and the condensate. However, this 
process is sensitive to the momentum anisotropy of the 
gap, and its rate .vanishes for an ideally isotropic gap, as 
noticed long agot^l in the context of the branch imbal- 
ance relaxation in NS junctions. On the contrary, in a 
d-wave superconductor, such a charge relaxation occurs 
at a time scale of order the impurity scattering time. As 
a result, compared with a normal metal (or an s-wave 
superconductor) , the quasiparticle charge diffusion mode 
is missing in a d-wave superconductor, which reduces the 
constant C from four to three. This reduction is a robust 
many-body effect of the anisotropic pairing symmetry. 

The plan of the paper is as follows. First, I show that, 
in an s-wave superconducting grain, the level correlations 
are essentially the same as in the normal state. Then 
I show how the gap anisotropy eliminates quasiparticle 
charge diffusion in a d-wave superconductor, and outline 
the corresponding microscopic calculation. Finally, I de- 
scribe the applicability range of the result and its relation 
to the previous findings, and discuss the possible further 
developments. 

S-wave superconductor. Consider an s-wave super- 
conducting grain. In the approximation of a spatially 
uniform gap, superconductivity can be described as pair- 
ing of the exact time reversed eigenstates of the under- 
lying metal.Ej Thus, the exact quasiparticle energies E n 
in the superconducting state can be expressed via the 
exact quasiparticle energies e n in the normal state as 
per E n = y / e 2 l + A 2 . Therefore, the exact density of 
states vs(E) = J2 n S[E — y / e 2 +A 2 ] in the supercon- 
ducting state is simply related to the exact density of 



f 



states vn{s) = J2 n ~~ e «] m ^ e norm al state: 



M E ) = 



E 



y/E 2 - A 2 
\/£ 2 - A 2 



^[v^ 2 - A 2 - e„] 



(V^ 2 - A 2 ). 



As a result, the density of states (DoS) dimensionless 
autocorrelation function R$(E,E') = { ^e)Hus(e% - 1 
in the superconducting state can be expressed through 
the DoS autocorrelation function in the normal state 
Rff (e, e') = y^T^^l - 1 in the form 

•5 V J I (u N (e))(v N {e')) 



R$(E,E') = R?(y/E 2 - A 2 , V(£') 2 - A 2 ). 

From this simple argument, it follows that the level 
correlations in an s-wave superconductor are identical (up 
to the change of variables) to those in the underlying 
normal metal and that, therefore, the diffusion modes in 
the two systems are the same. 

Before moving further, it is instructive to classify the 
quasiparticle diffusion modes more precisely. As a two- 
particle process, diffusion amounts to a coherent prop- 
agation of a particle and a hole, carrying spin-^ each. 
These two spins can add to form a singlet, which cor- 
responds to the charge degree of freedom - or a triplet, 
which corresponds to the three spin degrees of freedom. 
Both in s- and in d-wave superconductor, the condensate 
carries no spin, and thus a quasiparticle cannot exchange 
spin with the condensate. As a result, in a disordered 
spin-singlet superconductor, the quasiparticle spin dopa 
propagate diffusively, as it was recently re-emphasized .Ell 
By contrast, the situation with the quasiparticle charge is 
more delicate. Perhaps the simplest way to observe this 
difference is by inspecting the Bogolyubov quasiparticle 
creation operator: 



7p T = UpC+^+VpC^, u 2 p 
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Impurity scattering is elastic, i.e. it conserves the quasi- 



particle energy E p — y e 2 + A 2 . In an s-wave super- 
conductor with the uniform gap, A p is a constant and, 
in the absence of the Andreev scattering that turns e p 
into — s p , the energy conservation implies conservation 
of Up and v p . This means conservation of the particle- 
hole content of a quasiparticle and leads to the effective 
charge conservation - even though a Bogolyubov quasi- 
particle, being a superposition of a particle and a hole, 
does not have a well defined charge quantum number. 
The same conclusion can be reached by considering di- 
rectly the expectation value of the quasiparticle charge 
Qp ■ 

Qp = u 2 p (+l)+v 2 p (-l) = 



A 2 



In an s-wave superconductor, and in the absence of the 
Andreev processes, Q p is conserved by the impurity scat- 
tering, which leads to the charge diffusion pole. 

D-wave superconductor. By contrast, in a d-wave 
superconductor, the gap A p has strong momentum de- 
pendence. Hence, even in the absence of the Andreev 
scattering processes, Q p is not conserved by impurity 
scattering. In other words, in a d-wave superconduc- 
tor, elastic scattering does not conserve the moduli of 
the Bogolyubov factors u p and v p , thus changing the 
particle-hole content of a quasiparticle. Physically, this 
means that the impurity scattering leads to the exchange 
of charge between the quasiparticle subsystem and the 
condensate at the time scale of order the impurity scat- 
tering time - and thus to the absence of the charge dif- 
fusion pole. 

The same conclusion can be reached in a different and 
more formal way, by using the Ward identities, which 
are a consequence of the symmetries of the system.E3 In 
the Nambu notations, the BCS Hamiltonian of a d-wave 
superconductor reads 



H = / dr^ 



t~3e(p At 3 ) + r 3 e(j) — fJ-B + t 3 u 



dRdr&(R+ -)n A(E, r)*(ii - -). 



Here ^ = (ip^^i) is the Nambu spinor, Ti arc the 
Nambu matrices, R denotes the center of mass coordinate 
of a Cooper pair and r denotes the relative coordinate. 
The external fields are the vector potential A, the electric 
potential <fi, the Zeeman field B in the z direction, and 
u is the impurity potential. The pair field A(i?, r) has 
been chosen real for the sake of simplicity, and assumed 
to have the d-wave angular dependence. 

The Hamiltonian above respects the gauge symme- 
try \P — ► Ugty = exp [*73^Xg] * ( to b e accompanied 
by the appropriate change of the potentials and of the 
gap function), and the z axis spin rotations symme- 
try ^ — > U s ^f = exp [*^X»] * ( to be accompanied by 
the change of B)JtB Tha-.sought identities can be ob- 
tained in the usual wayJ23 by identifying the variation 
of the Green function under an infinitesimal transforma- 
tion G — > U GU + with the first-order perturbative cor- 
rection. In the absence of the external fields A, B and 
</>, and assuming spatially uniform gap (A(i?, r) = A(r)), 
the two symmetries imply two (Q = 0, fl = e — e' — > 0) 
Ward identities for the disorder-averaged Green functions 



G r,a(^p) 



e - - TiAp - t 3 s p , E fl (e) = S^(e) 



and the vertex renormalizations (tq)ra 
(A p T 2 )flA : 

- 2ir 3 S^(e) = (e - s'){t 3 ) ra + 2i(A p T 2 ) RA 



t 3 ) ra and 



RA- 



(2) 
(3) 



The Nambu matrix tq corresponds to the z-component 
of the quasiparticle spin, and the vertex renermaliza- 
tion (tq) R a describes its propagation. Similarly, t 3 is 
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FIG. 1. The ladder series for the renormalization of the 
vertices (rj). 



the quasiparticle charge operator, whose propagation is 
accounted for by (t^ra- The subscript RA means that 
the vertex joins a retarded and an advanced Green func- 
tion at energies e and e' .EH For the illustration purposes, 
the leading approximation for the vertex corrections (the 
ladder series) is shown in Fig. [3. 

Direct inspection of the ladder series shows that 
(A p T2)ra oc {tz)ra, which leads to the conclusion that 
the vertex (t$)ra remains finite as (e — e') — > 0. Thus, 
the quasiparticle charge does not propagate diffusively. 
At the same time, as seen from Eq. the vertex {tq)ra, 
corresponding to the z-component of the spin, acquires 
a diffusion pole: 



\Tq)ra 



i(£-s')' 



Notice that, in a metal, the charge diffusion mode re- 
appears, as seen by sending A p to zero in Eq. (||). 

The calculation. The microscopic calculation of the 
level correlations amounts to finding the DoS autocor- 
relation function K(e,e') = {v(e)v(e')) — {v(e)){v(e')) . 
The mean square deviation of the number of levels in an 
energy interval of width E is then given simply by 



(8N%) = [ dede'K(e,e'). 
Je 



The calculation of K(e,e') amounts to evaluating the 
Feynman diagram on Fig. as done inB, and has to be 
performed in the four-spinor Balian-Werthamer space. 
The technical difference with regards to a metal being 
that now the Green functions reside in the matrix space, 
whereas the impurity ladder resides in the space of direct 
products of the two matrices. The calculation for nodal 
quasiparticles in a d-wave superconductor leads to Eq. uL 
with C = 3, and its details will be published elsewhere.EBI 

The validity range. The applicability range of this 
calculation is set Jay the possibility to treat the gap as 
spatially uniform.E3 Hence the energy interval E should 
be much wider than the fluctuations of the gap. The 
latter scale is set by the perturbation of the gap due to the 
impurity potential. Neglecting the Coulomb interaction, 
the perturbation of the gap due to a single impurity has 
the formEJ 




FIG. 2. The Feynman diagram for the leading term in the 
DoS autocorrelation function. 



SA(R) ~ A [uN(0)} [AJV(O)] F(R), 

where A is the value of the gap far from the impurity, 
[uN(0)] is the dimensionless strength of the impurity po- 
tential, [A-ZV(O)] is the dimensionless BCS coupling con- 
stant, and F(R) is a function decaying to zero at the 
lcngthscale of order the coherence length. Scattering off 
5A(R) is the Andreev reflection off the inhomogencitics 

of the gap 

-l 



with the rate r . 



A/ £F ] 2 [AAT(0)] 2 « 

r~", where r _1 is the impurity scattering rate. Thus, the 
main result of this paper holds for E much greater than 
t7 , and for the levels separated from the Fermi energy 
by more than t^ 1 . At the same time, E must be much 
smaller than the Thouless energy E c = D/L 2 jB which 
requires t^ 1 <C D/L 2 
grain size by 



The latter inequality bounds the 



L < 



WW)' 



which is much greater than the coherence length. Note 
that, for nodal quasiparticles in a d-wave superconductor, 
both I and £ are functions of the quasiparticle energy e 
and scale as 1(e) ~ lA/e, which further increases the 
upper bound on the grain size. 

Before closing, it is instructive to put the main result of 
this work, Eq. (0) with C = 3, in context. Equation (|l|) 
(with different values of C) is commonly associated with 
the random matrix theory (RMT)Jij which furnishes 
very general and powerful phenomenological framework 
for treating jasriom systems, and allows a symmetry 
classificationliiMlJ of possible universality classes. How- 
ever, a crucial underlying assumption of the RMT is that 
all the matrix elements of the Hamiltonian (including the 
matrix elements of the gap) be random and drawn inde- 
pendently from a broad distribution. This requirement 
automatically rules out the possibility to distinguish su- 
perconductors of different pairing symmetry. 

By contrast, the present work studies the intermediate 
energy limit when the matrix elements of the gap may 
be treated as completely non-random, the randomness 
being restricted to the diagonal of the Bogolyubov-de 
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Gennes Hamiltonian. In this limit, the Hamiltonian is 
only a "partly random" matrix, with the matrix elements 
of the gap fixed by the pairing symmetry. As shown 
above, the level correlations in such a "partly random" 
ensemble are sensitive to the momentum anisotropy of 
the pairing and are qualitatively different e.g. for s- and 
d-wave superconducting grains. 

More importantly, this work shows that, in anisotropic 
superconductors, the condensate assumes the role of a 
"charge reservoir" coupled to the quasiparticle subsys- 
tem, and affects the level correlations. It would be inter- 
esting to study this problem further by explicitly includ- 
ing the dynamics of the condensate, especially in view of 
the cuprate superconductors as an obvious experimental 
object. 
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ence Foundation (DMR 91-20000) through the Science 
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